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GROUPS OF THE FUNDAMENTAL OPERATIONS OF ARITHMETIC 
By G. A. MILLER 


By adding n to z,, then this sum to x), ete., there results an infinite series 
of numbers*, n,n + 2, n+ 2x,,-- +. If these numbers are taken with 
respect to a given modulus (m) only m/d of them are distinct, d being the 
highest common factor of mand x,. These distinct numbers are of the form 
n+kd,k=1,2, +--+, m/d. If, represents the smallest one with respect 
to modulus m, their least positive residues are ng + kd, k = 1, 2, - - -, mid. 
These elementary principles of number theory will be frequently employed 
in the following developments. 

Instead of adding x to one number we may add it to a series of numbers 
Xj, %, +--+ , 2. If the sums thus obtained are again added to this series in 
any order, and these operations are repeated, the resulting numbers again be- 
long to m/d distinct classes with respect to modulus m, d being the highest 
common factor of 2,, %, - - +, 2, m. Their least positive residues are 
again given by ny + kd, kk =1,2, - - -, m/d, where ng is the smallest number 
with respect to modulus m that is obtained by the given operations. 

As the definitions of a group clearly apply to these operations and as all 
the numbers are obtained by a repetition of the same operation, it is conven- 
ient to say that adding to one or more given numbers and taking the results 
with respect to a fixed modulus gives rise to the cyclic group whose order is 
equal to the number of least positive residues of all the resulting numbers. For 
instance, when n = 8, x, = 15, m = 20, the least positive residues of the resulting 
numbers are 3, 8, 13, 18. These become the cyclic group of order 4 when 
they are increased successively by some multiple of d=5. These results 


* Allthe numbers employed will be regarded as positive integers unless the contrary is 


stated. 
(41) 
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remain unchanged by adding x, = 35 to the given conditions, but if x, = 12 
the cyclic group becomes of order 20 since d=1 in this case. 

These considerations are so elementary that they offer little interest in 
themselves. They are explicitly stated, partly for the sake of completeness 
and partly on account of their fundamental importance in the following devel- 
opments. It should be carefully observed that we add only to 2, %, - - -, 
zx, If the resulting numbers were added together, the order of the cyclic 
group would be m/d’, d' being the highest common factor of d andn. The 
fundamental thought in what precedes is that the groups of addition are cyclic. 
We proceed to consider subtraction. It will appear that its groups are of the 
dihedral rotation type. 


1. Groups of subtraction. If we subtract x from 2, the resulting 
number (x, — n) may be supposed to take the place of either x, or x. In the 
former case subtraction is an operation of infinite period unless we consider 
a modulus; in the latter case it is of period two, since x, — (%,—”) =n. 
We shall here adopt the latter point of view and hence regard subtraction 
from x, as a reflection on the point x,/2, the numbers being represented in the 
usual manner by points on a line. 

Subtracting x from x, and then subtracting the resulting number from 2, 
is equivalent to reflecting n on the point x,/2 and then reflecting x, — n on the 
point x,/2. The result is equal to a translation through x, — x,; that is, n 
becomes n + (%,—2,). It was seen above that these translations with respect 
to modulus m generate the cyclic group of order m/d, d being the highest com- 
mon factor of x,— 2%, and m. Hence the operations of subtracting from x, and 
x, generate the dihedral rotation group of order 2m/d.* One of the most 
interesting examples along this line is furnished by taking the complement 
and the supplement of an angle. As these two operations are equivalent to a 
rotation through 7/2 they generate the group of order 8 which is identical with 
the group of movements of the square.t It contains two operators of period 4, 
five of period 2, and the identity. 

While the general theory of the group obtained by subtracting from two 
numbers follows directly from the theory of the groups which are generated 











* Bulletin of the American Mathematical Society, vol. 7 (1901), p. 424. 
J Davis, Nebraska University Studies, vol. 4 (1904), p. 231. 
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by two operators of order two, yet it seems desirable to add a few details 
with respect to the special problem on hand. A reflection followed by a 
translation is equivalent to a reflection of the original point on the point 
obtained by moving the original point of reflection through half the transla- 
tion. Hence a reflection followed by a translation is always of period two, 
and a reflection followed by a translation followed by the same reflection is 
the inverse of the translation. In other words a translation is transformed 
into its inverse by any reflection. Hence the order of the group generated by 
any translation and a reflection is twice the order of the group generated 
by the translation alone. When the translation is equal to two given reflec- 
tions the group generated by the translation and one of these reflections 
includes the other reflection. 

When n is subtracted from a series of numbers x,, %, - - - , x! and the 
remainders are again subtracted from these numbers, in any order, these 
operations of order 2 generate the dihedral rotation group of order 2m/d, 
where d is the highest common factor of m, and of all the differences of the form 
x, — 2g. This follows immediately from the fact that the translations which 
are equivalent to the pairs of these reflections generate the cyclic group of 
order m/d. As there are just two positive points less than m with respect to 
which all reflections are equivalent modulo m, there are just 2m/d different 
values of n modulo m which go into less than 2m/d distinct numbers by the 
operations of the group of subtraction. Each of these special values goes into 
m/d distinct values, differing by multiples of d, under this group: 7. e. these 
2m/d special values constitute two sets such that all of those of one set go into 
each other under the group.* Congruent numbers modulo m go into con- 
gruent numbers under all of these operations. 


2. Groups of multiplication and groups of division. Since 
the groups of multiplication have received considerable attention f it seems 
necessary to do little more here than state those results which are funda- 
mental in the study of the groups of division. There is, however, one impor- 


* For example, by taking the complement and the supplement we obtain 8 distinct angles 
from a given angle provided it is nota multiple of r/4. In the latter case only 4 distinct angles, 
differing by +/2, are obtained hy these operations. 

+ Cf. ANNALS OF MATHEMATICS, ser. 2, vol. 2 (1901), p. 77. 
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tant principle connected with these groups which has not received sufticient 
attention: viz. the necessary condition that a set of numbers becomes a group 
when they are combined according to multiplication. It has been observed 
that a sufficient condition is the fact that the set is composed of all the num- 
bers which are less than a given modulus (m) and prime to m. That this 
condition is not necessary follows directly from the fact that the four numbers 
2, 4, 6, 8 become a group when their products are taken with respect to 
modulus 10. 

The necessary and sufficient condition that a given set of numbers, which 
are distinct with respect to modulus m, generates a group with respect to 
multiplication is that the highest common Juctor of the modulus and a number 
of the set is independent of the choice of this number, and that m divided by 
this highest common factor is prime to this factor. That this condition is 
sufficient follows directly from the fact that by dividing by the highest com- 
mon factor we obtain a set of numbers each of which is prime to the new 
modulus. The eondition is also necessary to insure that every number of the 
set can be obtained by combining any other number of the set with numbers 
of the set. For example, 5, 35 generate the four-group 5, 25, 35, 55 mod 60; 
while 5, 25 generate the group of order 2 with respect to the same modulus. 
In both of these groups 25 corresponds to the identity. In general, the 
number which is congruent to unity, modulo m divided by the common 
factor, corresponds to the identity. 

When the elements of a group are not prime to the modulus the group 
is identical with the one obtained when these elements and the modulus are 
divided by their highest common factor. In other words, all multiplication 
groups can be represented by elements which are prime to the modulus. 
Hence this special class of abelian groups is identical with the groups and the 


sub-groups of the groups of isomorphisms of cyclic groups, whose properties 
may be assumed to be known.* We proceed to the groups of division and 
shall find that these are obtained by simply adding operators of period two to 
the groups of multiplication. 

The period of the operation of division as that of subtraction depends upon 
the point of view. In each case the group may be regarded as the same as 








*ANNALS OF MATHEMATICS, loc. cit. ; Lucas, Théorie des nombres, 1891, p. 395. 
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that of the corresponding direct operation. That is, we may regard subtrac- 
tion as an operation on the minuend so that the remainders become new 
minuends. From this standpoint there is no difference between the group 
of subtraction and that of addition. Similarly, if we regard the quotients as 
new dividends there is no essential difference between the group of division 
and that of multiplication. If, however, the quotient is regarded as a new 
divisor, the operation of division is always of period two, instead of being of a 
finite period only when a modulus is employed. In what follows division will 
be regarded as of period two. 

Geometrically the operation of dividing x, by n may be regarded as an 
inversion of n followed by a stretching in the ratio x, That is, any point » 
goes into the point x,/n by this operation. * It is geometrically evident that 
this operation is of period two. By dividing two distinct numbers x,, x, there 
result two operations of period two whose group will now be considered. 
The product of these two operations, in order, brings n into (x/x,)n. If each 
of the numbers 2, x, is prime to m there will be just one number (;) less than 
m such that x2, = x2, 2; mod m.f We shall regard x,/x, and x, as equal with 
respect to modulus m. Hence the powers of x,/x, constitute a cyclic group 
modulo m and the two operations x,/n, x,/n generate a dihedral rotation group 
whose order is twice the order of this cyclic group. 

The general case of dividing the / distinct numbers 2, x, - - - , 2 by 
presents scarcely any new considerations. If m is prime to all the numbers 
Ly, Lg, + + + 4 9, the group (G) generated by these/ operations of period two 
is obtained by adding operations of order two to the group (//) of multipli- 
cation generated by the quotients 2,/x,, 2/73, - - - , %2/%, - + + 5 %4/%. 
The order of @ is twice the order of //, and the additional operations of order 
two transform each operator of // into its inverse. That is, this group 
consists of an abelian group, which is found in the group of order ¢(m) gen- 
erated by all the numbers less than m and prime to m, and as many more 
operators of period two as there are in this abelian group. 

If one of the numbers 2, 2%, - - - , 2, occurs among the quotients 2x/x,, 








* This clearly remains trae even when the numbers are complex. 

+ This is also true when m, zx; and m, x2 have the same highest common factor (d) and 
m/d is prime to d, as is evident from the preceding developments. For brevity we shall not 
mention this evident extension in what follows. 
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1/2, + + + 5 %y/®yy- + + 5 1/2, then all of them are found among these quo- 
tients. Hence the coeflicients of n and 1/n involve either g or g/2 distinct 
numbers modulo m, g being the order of (, and n goes into g or g/2 distinct 
values under G as x, does not or does belong to //. For instance, by divid- 
ing 2 and 6 by nmodulo 11 the following ten distinct coetlicients result :* n, 4n, 
2 2 2 2 2 

dn, 9n, 3n,- + £9 es ee 
er ’ 6 2 6 7 2 ee ee 
ing 6 by these numbers since tot tas etc. For special values 
of n these numbers become a group modulo m only when x is equal to one 
of the given coefficients. The order of this group is equal to the number 
of distinct coefficients. The prominent thought of this section is that the 
group of division is obtained from the corresponding group of multiplica- 
tion in the same way as the group of subtraction is derived from the corre- 


No new numbers are obtained by divid- 


sponding group of addition. 


3. Groups of subtraction and division. Since the operations of 
dividing x, by n and of subtracting n from x, (2,/n, 2; — 7) are of period two, 
it is of interest to inquire into the group generated by them. As is well 
known the six values of the anharmonic ratio are obtained from any one of 
them by a repetition of these two operations when x, = 1.f Similarly, each 
of the six trigonometric functions (sin?x, cos?«, — tan?x, — cot? xz, sec? x, 
ese? x) is obtained from any one of them by the same operations. That is, the 
operations of inversion and subtraction from unity generate the symmetric 
group of order six. Since these two operations performed in succession are 
equivalent to an operation of period three, this result follows directly from 
the theory of the dihedral rotation groups. 

When x, = 2, the given operations lead to another well known group ; 
viz., the group of movements of the square. This follows directly from the 
fact that the two successive operations 2/n and 2 — n are equivalent to the 


; 2 a p 2 
operation 5 which is of period 4. In fact, 5) = 2/(2 - >=) 
2—n 





—Nn 2—n 








S .. 2—n\? 
i <d and ( ~) =n. When 2, = 3, there results the dihedral group 
—n l—n 





* For convenience the unknown (n) is also given. 
tCf. Burnside, Theory of groups of finite order, 1897, p. 18. 
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, Poy : . . 
of order 12 since 3, 8 of period six.* It is easy to prove that this group 


is of infinite order for every larger value of z,. 
The last statement may be proved as follows: Let s, = x, — n, 8, = ~ 3 

n 

x, — (x%,—1)n 
l—n ; 
Ba, ~_ (a > aB)n _ a al B'n 
dx, — (y + dx,)n y+ On 





r 4 x q 
Then will s, 8. = 7, — a? (818)? = In general let (s,s,)* 


a+ Bn 


“a Then will (8, 8,)* +} = . If the value 





of each of the fractions 8/a, 5/ylies between — : and — 2 and if x, 2 4 then 


/ U 
will the fractions al’y? satisfy the same conditions. As these conditions 


are satistied in(s,s,)*, they must be satisfied for the higher powers. That is, 
the two operations in question generate a group of infinite order whenever 
x,24. For negative values of x, the ratios §’/8', &'/y' are negative and not 
zero whenever 8/6, 8/y satisfy these conditions. Hence s,s, is also of an 
infinite period for all negative values of x. 

From the same substitutions it follows also directly that s,s, is of an 
infinite period whenever z, is a proper fraction ; for, if a,8 involve the same 
highest power of x, and if the coefficient of this power is unity in each of 
these terms then will the same be true of a’, 8’. It is evident that these 
remarks apply equally to y,& and 7,8. It should be observed that the 
groups of this section are independent of any modulus while all of those which 
precede, with exception of the trivial group of order two, are dependent upon 
a modulus. 


* These examples are special cases of the following general theorems: whenever a= ty 
the operations of substracting n from x; and dividing x by n give rise to the symmetric group 
of order 6; when rn = 2x, these operations lead to the octic group, and when xj = 3zq they 








generate the dihedral rotation group of order 12. When x; = 0 while zy has any other given 
value, there results the group of division and change of sign. This is evidently the well known 
four-group. When z,, 2% do not satisfy any of these conditions, the two given operations 
cannot generate a finite group which transforms every rational point of the plane into a rational 
point. Hence these four-groups are of fundamental importance in the theory of point trans- 
formations of the plane. The examples given above illustrate these groups when 7 = Zp. 
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48 MILLER 


While it is well known that the same group frequently presents itself in 
many different investigations, yet it may be in place to emphasize the fact that 
the octic group is the group of taking the supplement and the complement of 
an angle, of dividing 2 by n and subtracting » from 2, of the movements oi 
solid space which transform a square into itself, of the three valued function 
in four letters, of the substitutions in four letters which transform the cyclic 
group in these letters into itself, of isomorphisms of itself, of subtracting n 
from 3 and 6 modulo 12, of dividing 1 and 2 by n modulo 5, ete. The num- 
ber of such illustrations is always infinite but in this case some seem to be ot 
especial interest. . 








LINEAR DIFFERENTIAL EQUATIONS WITH DISCONTINUOUS 
COEFFICIENTS 


By Maxime Bocuer 


Ir we consider the problem of determining the simple harmonic vibrations 
of a stretched wire whose cross section changes abruptly at some points, or 
the analogous problem of the flow of heat in a rod of similar character, we are 
led to a linear differential equation of the second order whose coeflicients are 
discontinuous functions of the real variable «. Such “finite discontinuities,” 
or even much more complicated ones,* in the coefficients of the differential 
equation should not be regarded as being in any true sense singular points of 
the equation, the only peculiarity which they produce in the solutions being a 
discontinuity in the second derivative. The same is true when the coefficients 
of the differential equation have discontinuities where they fail to remain finite, 
provided merely that the integrals of the absolute values of the coefficients 
converge when extended up to such points of discontinuity. The proofs of 
the fundamental theorems concerning linear differential equations can be given 
in such a form as to cover all these cases, and, a fact to which I wish espe- 
cially to call attention, this greater generality does not involve any material in- 
crease in complexity or difficulty of the proof. In the following pages I have 
given such a treatment, contining myself for the sake of simplicity to equa- 
tions of the second order though the method used admits of immediate and 
obvious extension to linear differential equations of any order. 


1. General preliminaries. We shall be concerned exclusively with 
real values of 2 which lie in the interval ° 
(7) asa2csh. 

While some of the functions with which we shall deal are continuous 
throughout (7), others are of a more general character. These we designate 
as integrable according to the 











* As for instance the discontinuity which sin 1/2 has at the point z =0 
(49) 
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Derinition. A function f(x) ts said to be integrable throughout (1) if it 
is continuous except at a finite number of potuts, and of Jy x)ydx converges 


when extended over any portion of ( T).* 

The following facts concerning integrable functions are well known, or 
‘an be easily proved : 

1) If d(x) has only a finite number of discontinuities and (2) is in- 
tegrable, then (2) is integrable. If this is the case we shall speak of ¢ as 
absolutely integrable. 

2) If od). dy, -- + d, are integrable (or absolutely integrable), the func- 
tion ¢, + d, ---- - d& is integrable (or absolutely integrable). 


3) If d(.) is integrable, | o(x)dx is continuous. 


4) If d(.r) is absolutely integrable and f(z) is continuous, /-¢ is abso- 
lutely integrable. 

+) If b() is integrable and f(2) is continuous and has a continuous (or 
absolutely integrable) derivative, 7-6 is integrable. (Proof by integration by 
parts. ) 

It is in no way necessary, except in §5, to require that the functions with 
which we deal be real, although, as has been said, the independent variable x 
is supposed real. By slight and obvious changes the whole treatment can be 
carried over to the case of analytic functions of complex variables. This has been 
indicated briefly in a simple case in $4, where one of the independent variables 
(namely 2) is allowed to become complex. Only slight changes would be nec- 
essary in order to allow « also to become complex. 


2. The Existence Theorem.+ We will consider the ditferential 
equation 


dey : dy 
(1) dg TPO) a + Ue) Y = (4), 


*It should be carefully noticed that this definition is more general than Riemann’s in thag 
it does not require the function to remain finite in (J), but less general in that it admits only a 
finite number of points of discontinuity. . 

TC]. the earlier papers of the writer: Bulletin of the American Mathematical Society, March, 
1899; American sournal of Mathematics, vol, 24, 1902, where this same question is treated by 
slightly different methods, 
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and assume that r is an integrable function throughout (7) and p and q are 
absolutely integrable there.* 

The question now arises what we should understand by a solution of (1), 
the difficulty being that we can clearly not demand that a function which satis- 
fies (1) at other points should also satisfy it when one or more of the fune- 
tions p, g, 7 become discontinuous. If we make no demand whatever on the 
solution at such a point of discontinuity of the coefficients of (1), there will 
clearly be no connection between the values of the solution on the two sides 
of this point. It is desirable for most purposes to demand that the solution 
itself and its first derivative be continuous at such points. Thus we lav down 
the following 

DEFINITION. By a solution of (1) ts wnderstood a Junction of « which 
at every potnt of (1) is continuous and has « continuous first derivative, and 
at every point of (1) where p,q, rare continuous has a second derivative and 
satisfies (1 a 

We now establish the existence theorem : 

If cas any pomt of (/) and Ys ny coustants, there exists a solution 
of (1) which satisfies the ausiliary conditions 


(2) (ec) = ¥, y'(e) = ¥)- 


Let us establish this theorem first in the special case in which p and q¢ 
vanish at every point of (7). The general solution of (1) is then 


y =| ( renae) ae { Ar —c)+ Bb, 


a function which is obviously continuous and has a continuous first derivative 
throughout (7). If we let l= y,, B= x, conditions (2) are satistied. Thus 
the theorem is proved in this special case. The expression for y which we 
have obtained can be thrown into a more convenient form by an integration 


by parts, namely : 


Y= | r(x)yde - | wr (ayde +y(%—e) + y¥- 


*It is not unlikely that it would be sufficient to require that p and q be merely integrable 
The writer would be glad to learn of « method by which this more general case could be 


treated. 
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By a slight change of notation we get the following expression alongside of 
which we write for reference the value of 4’: 


ly =['e ~ £)r( dE + n(2 -— 4% 
(3) : 
\y=| r(E\dE +. 


We now treat the general case by the method of successive approximations. 
Starting from the value y, = 0, we compute the successive approximations 
Yi, Yo. +--+ bY means of the recurrent relation : 

Yn = — PYn—1» — Wars ‘e (n=1, 2,--+)y 
and the auxiliary conditions 
- ' 
(CO) = ¥ Yn) = Ne 


We shall prove that 7, approaches a limit as n becomes infinite, and that this 
limit is the solution of (1) whose existence we wish to establish. 
If we write 
u(x) — Yn > (xz) — Yrl“)s 
¥, is obviously the sum of the first n terms of the series 
(4) Ug(4) + U(r) + Wy(2) 4+ -- 


We must then prove that thix series converges and represents a solution of 
(1) which satisties conditions (2). It is convenient to consider by the side 
of (4) the series 

(5) U(L) + Wy(e) + UZ(L) + -- 


Since uv, and y, are identical, we see that up and uj are given directly by 
formule (3). The subsequent u’s may be determined by the recurrent relation 


u, = — p(*) Wer — YL) ny (r= 1, 2, vee) 
with the auxiliary conditions 
u,(c)=0, uw (c)=90, 


as we see by subtracting from one another the recurrent relations and also 
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the auxiliary conditions for determining y, ,, and y,. Applying (3) we now 
tind 

gq (M@=-[ © -PPOma® + 1) a, 
(6) “ 
une) = =f (PE h(E) + 118) nan (B)) AE 


Thus we see that all the terms of (4) and (5) are continuous functions of x 
throughout (/). 
Let us now introduce a positive constant / subject to the conditions 


(7) tal, Iz hea, 


We then infer from (6) the inequality 


(%) my f SEL (ine) +10!) (ie) | + ea BD!) |e 
(n= 1, 2, ---). 
Since uy(x) and uo(x) are continuous throughout (/), there exists a posi- 
tive constant C’ such that 
(9) u(x) = C, ug(x) = 


We are now in a position to prove the fundamental inequalities 


— I B [(2® — ) dt | 
nm. 


u(x) 
Since (10) reduces to () when n = 9, it is sutticient to prove that if it is true 
when n = A, it is true when n = & + 1. 

Letting » =’ + 1 in (8), and replacing |(u,(&) and |u,'(&)| in the 
second member of (8) by values obtained by letting n = * in (10), we get 


py. (2) 


Uy. +1 (x) 


so Of (re) + 9(€) LL (p<) +/9(&) ) dé ] dg|. 


From this we deduce the special case of (10) which we wish to establish 
by making use of the formula 


food (dé) dé= milf “$Edé) 


whose correctness follows by an integration oh parts. 





(10) 
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If now we let 


M= [ (re + 4(&) ) dé. 


“ 


we obtain from (10) the simpler inequalities 


u(x) ) 


(11) w(x) | | 


These inequalities show that at every point of (/) the terms of the serie 
(4) and (5) do not exceed in absolute value the terms of the series 


(' (20.0)? i 


»' 


and, since this is a convergent series of positive constant terms, it follows by 
Weierstrass’s fundamental test that the series (4) and (5) are uniformly con- 
vergent throughout (/), and therefore, since their terms are continuous, rep- 
resent continuous functions throughout this interval. Moreover, by a well 
known theorem concerning uniform convergence, we see that (5) represents 
the derivative of (4). 

Denoting by y(«) the function represented by (4), it is clear that ¥ 
satisfies conditions (2). In order to prove that it is the solution of (1) whose 
existence we wish to establish it is then merely necessary to prove that 7 satis- 
fies (1) at every point of (7) where p, 4, 7 are continuous. 

Let x) be any such point, and surround it by an interval (-/) lying within 
(7) throughout which p, 7, r are continuous. If we multiply (5) by — p(x), 
(4) by — 7(«) and add the two resulting series together, we get, as we see by 
referring to the differential equation for w,,, 


(12) — p(%) y'(x) — 4( 2) y(r) = uy" (2) + y"(2) 


From the way in which this series was obtained it is clear that it is uniformly 
convergent throughout (-/) and that its terms are continuous there. Accord- 
ingly, since it may also be obtained by differentiating (5) term by term (the 
first term being omitted) its value is 


y'(£) — ug (a2) = y"(v) — r(x). 


Substituting this value in (12), we see that y satisties (1) throughout (./). 
That the solution whose existence we have established is unique is stated 
by the following theorem which we now proceed to prove : 
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If two solutions of (1) satisfy conditions (2), they are identically equal 
throughout (/). 

The difference of two such solutions is a solution of the homogeneous 
equation 


d 7] dy 


1s —og = r)— ola = 
( ) Fe pt as yiejyy =, 


( 
satisfying the auxiliary conditions 
(14) y(e) = 9, yi(c) = 0. 
Hence our theorem will be established if we can prove the following more 
special result : 
A solution y, of (15), which satisfies the auxiliary conditions (14) vanishes 
at every point of (1). 
We prove this by a method due to Sturm.* Let #9 be any point of (7), 
and consider a solution y, of (13) satisfying the auxiliary conditions 
Y2(o) =O, ¥'x(9'9) = 1. 
By Abel's formula 
' . : a pdx 
Wor) = n(")ye (2) — w()y (24) = he os 
where 4 is a constant. Since Wc) = 0, we see that = 0. Accordingly 
Wsro) = 0. But Wao) = ¥, (7%). Thus we see that 7; vanishes at 2» which 
was any point of (/). 


3. Equations which involve parameters.; It often happens that 
the coeflicients p, 7, r of (1) involve besides the variable x one or more parame- 
ters, and it is then important to know the nature of the dependence of the solu- 
tions of (1) on these parameters. In what follows we confine our attention to the 
case of a single parameter, but the methods and results can be immediately ex- 
tended to the case of any number of parameters. 

We consider the differential equation 


d *y 


di 


(15) 


* Liouville’s Journal, vol, 1 (1836), p. 109. 
i Instead of allowing the coefficients of (1) to vary through the introduction of one or more 
parameters, we may, without any essential complication, adopt the much more general stand- 
point of the calculus of variations. C/. the paper in the American Journal already referred to. 
The conditions imposed there, at the bottom of p. 315, on the variation of the coefficients might 
easily be made much less restrictive, as will be seen by comparison with the present paper. 


ee ei ee ee 


ee 
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where the real parameter is confined to the interval 
(L) 1, =r Shp. 


We shall speak of the pair of values (x,’) as a point of the region (J, L) 
it x is a point of (7) and Xa point of (1). We impose on the coefficients of 
(15) the following conditions : 

(A) For every value of X in (L) r is integrable and p and q are absolutely 
integrable throughout (I). Furthermore the function 


[ rx ryaz 


is a continuous function of the two independent variables (x, X) at every point 
of (I, L). Finally, \ being any fixed point and r, 4 variable point in (L): 


ie 
| p(#, %) — p(2, Ay) dx=0, 


16 
aad, lim 


*b 
Ag = A G(s Ax) — 7(4, Ay) dc = .* 


The fundamental theorem is the following where we append an accent to 
a function of two variables to denote differentiation with regard to the first 
argument. 

Tf y(X) and y,(A) are continuous functions of X throughout (L), and ¢ is 
any point of (I), and the coefficients of (15) satisfy conditions (A), then the 
solution y(x, ) of (15) determined by the auxiliary conditions 


y(e,r&) =A), y'(e,&) =NOA) - 


is a continuous function of the wo independent variahles (x, ) throughout 
(I, L); and the same is true of y'(a7, X). 

This theorem will be proved if we can show that all the terms of the 
series (4) and (4) are continuous functions of (x, ) throughout (/, 1), and 
that these series converge uniformly in this region. 








*It can readily be shown that when conditions (16) are fulfilled the integrals 


S’ px, d)dz, S'q(x, dz 
are continuous functions of (+, \). It is an open question whether (16) might not be replaced 


by the less restrictive conditions that the two integrals last written be continuous functions of 
(4,2). Cf. foot-note p. 51. 
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In establishing the continuity of these terms it will evidently be sufficient 
to consider the functions u{ since the functions u; may be written 


Up = fw dz + y(r), 


Zz 
, . ‘ 
ed (é = I, 3, +++), 
JC 


so that their continuity as functions of (7, X) follows from the continuity of 
the ws. 
Since the function 1 is given by the second formula (3), it is obviously 
a continuous function of (2, %) on account of the restriction placed on 7 in 
condition (A). Accordingly we use the method of mathematical induction, 
assuming that it has been proved that ”,_, and w!,_, are continuous fune- 
tions of (x, 4) and then proving that wv! is continuous. For this purpose we 
use the value of u/, given in the second formula (6) which may be written as 
the negative of the sum of the two integrals 


(17) [re Aju’, —1(€, r) dé, [ag A)u,_1(&, A) dE. 


We wisb to prove that these two integrals are continuous functions of 
(x, X). Denoting the first of them by F’(7,A) we must then show that 


AF = F(x + Ax, ’ + Ad) — F(x, A) 
approaches zero as its limit as Az and Ad approach zero independently of each 
other. Now we have: 


“z+ Az 
aF =| [ p(—, X+ Ar) — pc&, A) Ju’,_1(E, A + AAAE 


“r+ ar 
+| pcé, A) [u', 1 (€. A + AA) — 1 (&, A) €, 


fz + az 
+| pe, Aju’, 1 (&, r) dé. 


We will denote these three integrals by ¢), ¢:, $3 respectively, and we 
will show that they all approach zero as Ax and Ad approach zero. 
Since u/,_, is continuous, a positive constant Iv can be chosen so that 


throughout (/, L) 


| uy 4 (Z, r) = i. 






ee eee 
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Accordingly we have 


$/<K | "| p(é, A+ Ad) — p(Es 2) | [ae 
» 
=f | p(& A+ AA) — pr), dé, 


and this is seen to approach zero by (16). Furthermore, 


‘ “b 
>: | P(E A) [un a(E, A + AA) — wa (E, A) | dE S | | p(&, ») | dé, 


where » is the greatest value of «},_\(,% + AA) — u),_,(#, A) in (7). 
Since » approaches the limit zero as AA approaches zero, ¢, does the same. 

Finally in considering ¢; let us notice that since p is absolutely integrable 
throughout (7) and w/,_, is continuous there, pu},_, is absolutely integrable 
there, and hence its integral ¢, is a continuous function of the upper limit of 
integration, ¢. e. of Az. This shows that ¢, approaches the limit zero with Az. 

We have thus proved that the first of the integrals (17) is a continuous 
function of (2,). The same proof applies to the second if we interchange p 
and g, and w)_, and u,_,. The proof that all the terms of (4) and (5) are 
continuous functions of (7, A) is thus complete. 

To show that these series converge uniformly in (J, L) we have merely 
to repeat with very slight changes the reasoning given in §2 for establishing 
their uniform convergence in (£). Thus the constant C’ must be so chosen 
that (9) holds throughout (/, L), and this is possible since uy and w', are now 
continuous functions of (2, X) throughout this region. Furthermore 


| ( P(E, AY) + 1G A) ) de 


is now a function of X. It is easily seen by (16) to be continuous throughout 
(ZL), and therefore we can introduce a positive constant .W so that 


M 2 ["( plé, A) + | 9(& A) jae. 


Formula (11) will then hold and the proof of uniform convergence is complete. 

The condition (1) which we have imposed on the coefficients P> Y 7 is so 
general that it is worth while for us to notice some simple cases which 
frequently occur in practice and which are included under it. We mention 
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here three such cases each of which is included as a special case under the next 
following one: 

(a) Condition (A) is fulfilled if p, g, r are continuous functions of (x,r) 
throughout (J, 1). 

(4) Condition (A) is fulfilled if », 7, rare continuous functions of (x, 2) 
throughout (/, 1) except when # = ¢, ¢, ---¢,: and a positive constant M 
exists such that except at the last mentioned points |p <M, \qg| <M, |r\< M. 


(c) Condition (A) is fulfilled if except when x = ¢, c, +++ ¢% 


p =S(®) $)(%,r), Y =S(©) $:(*r), r=f (7) $3(x, r), 


where $,, $2, $; are continuous functions of (7A) throughout (7, ZL); and 


I(x) = (# — |)™(H — G)™ +++ (@— )%, 
a), - ++ a being real constants greater than — 1. 


4. Equations whose coefficients are polynomials in \. By 
imposing further restrictions on p, 7, 7. y. y% We could make the functions y 
and y' have certain further properties, for instance that of possessing first 
partial derivatives with regard to X which are continuous functions of (7, d). 
Without stopping to treat these intermediate possibilities, we will consider what 
in @ sense is an extreme case, and yet one sufliciently general to cover the 
cases which most frequently occur in practice. The theorem is this : 

[fc is any point in (1), and y, y, any constants, and if p,q, r are poly- 
nomials in % whose coefficients are integrable (and in the case of: p and q ab- 
solutely integrable) functions of x throughout (1), then, y(x.r) being the 
solution of (15), which satisfies the auxiliary conditions 

YaA=% WA ="N 
y(x, ) and y'(x, A) are continuous functions of (x, r) when x les i (J) 
und » has any value real or complex. Furthermore, «9 being any point of (1), 
y(xo, A) and y'(x, %) are integral analytic functions of r.* 

In order to prove this let us consider a region (A) in the complex A- 
plane consisting of all the points within and on the boundary of an arbitrarily 


chosen circle. The functions p, 7, 7 obviously satisfy condition (A) of §3 if 


in this condition the region (A) is substituted for the interval (Z). Accord- 
ingly we see, from the proof there given, that (#, X) and y'(.r, ) are contin- 


* That is, analytic functions having no finite singular points. 


a 


oe te rr er nen ees 
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uous functions of (2, X) throughout (/, A). Moreover, the series (4) and 
(5) being uniformly convergent in (/, A), will, when we let x = 2, be uni- 
formly convergent in (A). From formule (3) and (6) we see that the terms 
of (4) and (5) are polynomials in %. We see therefore by a fundamental 
theorem of Weierstrass* that (a, %) and y'(#, X) are in (A) single valued 
analytic functions with no singular points. From these facts our theorem fol- 
lows since the circle (.A) may be taken as large as we please. 

Finally we note that the assumption made in the last theorem that 
and y, are constants is not essential. For if we denote by y, and y, the solu- 
tions of (15) which satisfy the auxiliary conditions : 


yi(e, A) = 1, y2(e, A) = 9, 
yy (ec, A) = 0, ¥3(C, r) 


the solution y which satisfies the conditions 


Yer ®) = AWA). -W'(Cy ®) = NCA) 
miv be written 
Y= VN + N Yr 


If then y and y, are single valued analytic functions withno singular 
points in a certain continuum (Ay) of the A-plane, we see at once that v(x, ) 
and y'(z, X) are continuous functions of (2, %) throughout (7, Ag), since by 
the last theorem this is true of 4), ¥{, Y2,. yj. For the same reason if x9 is any 
point in (7), v(x, %) and y'(.79, A) are single valued analytic functions with- 
out singular points in (Ay). 


5. Roots of functions which involve parameters. Let us con- 
sider a real function 6(x,) detined for all values of x in (7) when A belongs 
to a certain continuous or discontinuous range of values. We suppose that ¢ 
satisfies the following conditions : 

(a) @ has a first derivative with regard to x, 6'(7, X), which is continuous 
throughout (/) for all values of X with which we are concerned. 

(4b) As X approaches a certain value / (either continuously or in some 


*If in a connected two dimensional region of the complex plane a series converges uni- 
formly, and if throughout this region its terms are single valued analytic functions withou, 
singular points, then the same will be true of the function represented by the series. Cf. 
for instance Burkhardt, Analytische Functionen, 1st edition p, 139. 
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discontinuous manner as the case may be*), $(7, ) and ¢'(a, A) approach 
uniformly certain limits f(#) and f'(2). 

(c) The functions f(x) and /’(z) do not both vanish at any point of (/). 

(d) f(x) vanishes neither at @ nor at b. 

It should be noticed that it follows from (a) and (/) that J (x) and f'(x) are 
continuous throughout (7), and that f’ is the derivative of J. From (c) it 
follows that (7) can not vanish at an infinite number of points in (I), for 
otherwise these points would have at least one limiting point in (7) at which 
both f and /’ would vanish. 

TueoreM. J/f $(2x, X) satésfies conditions (a), (6), (ec), (ad) and if its 
limit f(x) has just k roots x, --+ x, in(1), then for all values of X sufficiently 
near tol, d(x, X%) has just k roots in (1), and these roots approach the values 
Zi, ° + + ®, as dX approaches 1. 


In order to prove this let us suppose the roots 7, - -- 7, to be arranged 
in order of magnitude (a <2, <2, <--- <2, <4), and consider the /: intervals 
(J) r—az,| £6 (¢=1,2,---k), 


where 6 is a positive real constant which we suppose taken so small that all the 
intervals (/;) lie in (/), and that no two of them have a point in common. 
Since f'(x) is continuous and, by (c), does not vanish at any of the points z,, 
5 may be taken so small that /’(2) does not vanish in any of the intervals (/;). 
This we suppose done. 

We define further the intervals -/ as follows : 


(Jy) asxszx,-— 5, 
(J) 4+6S8252,4,-—64 (¢=1,2,---AkA-—1), 
(-J,) x}. + 6 => 7 = h, 


Our theorem will be proved if we can show that no matter how small 6 
may be chosen, a positive ¢ can be taken so small that for every value of A (in 
the range with which we are concerned) for which 2» —/ < et, the function 
$(,A)has no roots in any of the intervals (-/;) and just one root in each of the 
intervals (J;). To prove this let us denote by m, the smallest value which 
J’(x) | has in any of the intervals(/;), and take € so small that when A —/) <€ 





*For instance we might have 7 = + © and the range of values of \ with which we are 
concerned might be the positive integers. 


tIfl = © we must replace this inequality by 1A: 


a 
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the inequality | ¢’(2, 4) — f(x) | < m, is satisfied. This makes it certain that 
$'(x, X) does not vanish in any of the intervals (/;) when |A—J| <e. Now 
denote by m; the smallest value of | f(x) in the intervals (J;), and take e still 
smaller than before, if necessary, so that when |A —/| < the inequality 
o(x,r) — f(a) |< my is fulfilled ; thus making it certain that, for these values 
of A,@(x, X) does not vanish in any of the intervals (-/;), and has the same 
sign throughout each of them that f(x) has there. But f(x) changes sign at 
every root z,, since f’(x;) + 0; accordingly, when A —/|< e, the sign of 
$(x, X) is alternately plus and minus in the intervals (Jo),(J;),---. Thus we 
see that $(a,), which has no roots in the intervals (-/;), has at least one root 
in each of the intervening intervals (/;). It cannot have more than one root in 
any interval (/;), for otherwise #’(x2) would vanish there by Rolle’s theorem. 
We will now apply the theorem which has just been proved to the solu- 
tions of linear differential equations of the second order, confining ourselves, 
however, to the case of homogeneous equations, since here the mere require- 
ment that the solution is not identically zero is suftictent to ensure its satisfy- 
ing condition (c). 
If in the equation 


Py 


dy 
+ p(x, r) + q4(%, A)y= 0 


p and q satisfy conditions (A) of §3, and if y(x, X) is the solution of this 
equation which satisfies the conditions 
Ya, %A)Y= WA), -W¥'(4, AY = NCA), 

where y and y, are continuous throughout ( L), and if, for all values of 
in (L), y(2, X) has at least kh roots in the interval a < x Sb which arranged 
a tehilin of smmgiltiills ee dnilite Sep 8), eka -- + Whom Hy - + my O00 
throughout (L) continuous functions of provided one or the other of the 
following conditions is fulfilled : 

(1) yA) does not vanish in (L); 

(2) (A) = 9, and y,(A) doex not vanish in (L).* 





*The object of these last two conditions is to prevent z,(\) from reaching the value @ 
This may be attained more generally by the following condition: 

(3) y and 7, do not both vanish at any point of (1), and the fraction 7,/y has, for values 
of \ for which 7 + 0, if any such exist, a finite lower limit, possibly negative. 

The cases mentioned in the text are, however, by far the most important. 
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We prove this theorem first in case (1). Here we have merely to prove 
that if / is any value of A in (L) the functions x,, - - - x, are continuous when 
X=/. Suppose that y(x,/) has just / + m roots in (J). If y(4,1) 40 
our theorem follows trom the one last proved which tells us that when 2 is 
sufficiently near to /, y(,) also has exactly 4: + m roots in (J), and that as 
X approaches 1, x,(A), - - - {4m_(A) approach respectively x,(/),- + - x4 ,(1). 
If, however, y(6, 7) = 0 we must extend the interval (7) a little beyond the 
point & by defining p and q for values of « a little greater than 6. For in- 
stance we may let p(x, 4) = g(z,A) = 0 whenz>b. If b'>b, y(',l) + 0 
since y(b, 1) = 0. Accordingly we may apply the first theorem of this section 
to the interval a = x = 8’, and we see that, for values of \ sufficiently near to 
1, y(a, X) has just & + m roots in this interval, and that these roots approach 
vy(1), +++ 24m (7) as their limits as % approaches /. Since / of these 2’s 
always lie in (J) by hypothesis our theorem is proved. 

It remains to prove the theorem when y(A)= 0. Here we determine a 
constant a’ > a but differing from a so little that y(x, X) does not vanish in 
the interval a < « Sa’ for any value of Xin (Z). This we do by taking a’ 
in such a way that y'(x, X) does not vanish in the interval a = x = a’, for then 
y which vanishes when x = a could not vanish again by Rolle’s theorem. By 
applying to the interval a’ = x = 4 the same reasoning which was before ap- 
plied to (J), the truth of our theorem follows. 
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SOME PHYSICAL SOLUTIONS OF THE GENERAL EQUATION 
OF THE V™ DEGREE 


By Rosnert E. Moritz 


SUCCESSFUL attempts to solve general equations by physical means date 
from about 1770. The various devices employed by Plato, Apollonius, 
Eratosthenes and others in the solution of the Delian problem did indeed afford 
solutions of the equation «° = 20°, but as is well known the Greeks never rose 
to a full conception of the general equation of even the third degree. 

In 1770 J. Rowning exhibited his Universal Constructor of Equations, 
constructed on perfectly general principles though the working model was 
limited to the solution of quadratics. During the century following, L. 
Lalanne (1840) succeeded in constructing algebraic balances for the solution of 
equations of as high a degree as the seventh. Since 1870 a great number of 
devices for the solution of the general equation of the 2" degree have been 
invented, based upon principles of mechanics, optics, hydrodynamics, elec- 
tricity and magnetism. A brief enumeration of these contrivances and their 
bibliography is given by R. Mehmke in the Encyhlopddie der Mathematischen 
Wissenschaften, vol. I, p. 1067. 

The object of the present paper is to present with clearness and in an 
elementary way several physical solutions which, it seems to the author, deserve 
to be more widely known, not only because of their intrinsic interest but es- 
pecially because of their pedagogic value. They offer a remarkable illustration 
of the inter-relation of the various elementary branches of mathematical science 
in revealing how various principles may lend themselves to the solution of the 
same problem. At the same time these solutions exhibit well the unifying 
effect of the mathematics, in that the most diverse physical principles are seen 
to conform to the same mathematical laws. Above all they form a good 
example of the correlation of all science, of how a comparatively difficult prob- 
lem may find its complete solution in apparently unrelated fields. These 

(64) 
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elements, together with the peculiar fascination of the subject, the simplicity 
of the physical laws involved and the elegance of the underlying mathematical 
principles, commend these solutions to a place in our elementary text books. 

While the ideas underlying the solutions as here presented are not new, 
the treatment is independent. The detailed description of the various devices 
has been omitted, for it is needless to say that the interest in these practical 
devices is essentially theoretical. Mechanical difficulties prevent the actual 
construction of working models for all but the lowest degrees. The best 
models as yet constructed give at best but a rough approximation to the roots 
sought. Besides it may well be questioned, whether since the perfection of 
approximation methods, such as Horner's, there exists any longer the need for 
actual solutions by physical means. For this reason the mathematical and 
physical principles of these devices have been emphasized rather than their 
mechanical execution. 

The first solution is Rowning’s,* based upon Segner’st method of evalu- 
ating polynomials. The second is due to Lillet and has been treated with 
modifications at some length by Cremona.§ The third solution is based upon 
Boys’|| as modified by Grant. The fourth, last in point of discovery, is 
Meslin’s**® hydrostatic solution. The fifth solution is that of Lucas,tt though 
the principle underlying it dates from Gauss. tf 


1. A kinematic method. It is evident that if a mechanism could be 
devised for. tracing the graph of the polynomial 


Yn = Agu" + ax"! 4+ ag"? 4+ --- +a," + a, [1] 


in rectangular coordinates, the distances from the origin to the points where 
this graph crosses the x-axis would give at once the real roots of the equation 


aya" + a,x"—! + agx"—* Se oe Ay an a, = (). [2] 


* London Transactions 1770. 

t Acad, Petropolis Novi Comment. 1758-59. 
t Comptes Rendus, 1867, p. 854. 

§ Graphical Statics, Oxford, 1890. 

|| Philosophic Magazine, 1886, p. 241. 

{ American Machinist, 1896, p. 824. 

** Comptes Rendus, 1900, p. 88S. 

tt Comptes Rendus, 1890, p. 965. 

tt Werke, vol. 3, p. 112. 
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4) a mechanism is made manifest from the following 


The possibility of suc 
method of constructing the value of , in [1] tor any given real value of x. 


FIG. 2. FIG. 3. 
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Draw the line ALY (Figure 1) and on it lay off AZ? equal to unity and 
AP equal to the given value of 2. Through A draw AY’ perpendicular 
to JLY, and through 7? and /? draw parallels to AY. 

On ut) take 


AA, =4,, A,A,_) = G@y—is °° *s Ay, = 4, A, A, = %, 


incasuring each distance in the positive or negative direction according as the 


sigan of the corresponding coeflicient is positive or negative. Through Ap 
avo line parallel to ALY, intersecting the parallels through P and 2 in 7, 

and /2, respectively. 

Draw the line 22.4, and let 2; mark its point of intersection with the 
line 177%). 

Draw 2/2; through 7, parallel to ALY intersecting 222?) in J?;. Draw 
yA, and let its point of intersection with 2, be P,. 

Similarly determine in succession the points 7, P, ---, P,—1, P,- 

The distance ?7?, represents the polynomial y,,. 

Figure 1 gives the construction for the case in which the coefficients are 
all positive and < 1. Figure 3 shows the construction when some of the coefti- 
cients are negative and w > 1. 
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The proof of the construction follows most readily from Figure 2, which 
is a repetition of Figure 1, with a more convenient notation. Through the 
several points A,, Ay, A;,---, A,_,, A,, draw lines parallel to AY, inter- 
secting /2F, in the point 2,, B,, B,,---, B,_., B, respectively. Moreover 
let» = Bho, yy, = Bk, +, y,-1 = BR, —1, Yn = RM, = PP, 


Now we need only to compare any pair of similar triangles such as those 


marked (7) and (//) in Figure 2 or 3 to see that for every such pair of 


triangles there subsists a relation such as 


Y2? 1=(y¥3 — 43): 2 
from which 
Ys = Yt + Ag. 


We thus obtain the set of equations 


Yo = lo, ) 
%; = Yor + 4, 
Ys = Ye + >, 
3 
> und 


Yuna = Yn—-2 + Anis 


Yn = Yui T any 





from which by successive elimination 
Y\ —_— yx i is 
Yq = Agk* + ayX + Ay, 


Ys = Ax? + aux? + ar + ay, 


Yu = Ax” + aya" 4+ age" 4 + + + Age + A), 
This last result may also be obtained otherwise by eliminating the y's in 
[3] simultaneously as follows. We multiply the first of the equations [3] by 
x, the second by x«"~', and similarly each equation whose first member is ¥, 
by 2"—*, and then add the results. The left member of each equation excey 
the last cancels against the first term of the right member of the equation jus: 
following and we have remaining 


Y= PP, = AyX" + gx" at ag*-* tees t+ Ay a. 
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From the foregoing construction a device for tracing the graph of [1] 
readily suggests itself. A), AN and RR, are replaced by fixed rods, of 
which A Y carries a graduated scale for the proper determination of the points 
A,, Ayia +++ Ay, Ap. BoP, Hi Pi, ete., oA, Ag, ete., and PPpo are 
replaced by adjustable rods properly grooved and pivoted so as to admit of 
parallel motion of the rod PP,. While this rod PP is made to move parallel 
to itself a pencil point at P,, traces a graph on a sheet of paper below. The 
points of intersection of this graph with the line AX furnish the required 
roots. 


2. A solution by visual means. A second solution is suggested 
by the following method of constructing the polynomial [1]. 
Let us first assume all the coefficients positive. 


Fi. 4. 




















Draw PA, = a (Figure 4), and from A, in a direction + 90° with PA, 
draw A, A; = Q). 

From <A, in a direction + 90° with A) A, draw A, Ay = ay. 

Similarly draw the remaining segments a5, ---, a, of the broken line 
PAA, --- A,. 

Beginning at P draw another rectilinear figure P,P, --- P,, having its 
vertices P,, Py, ---, P, respectively on the prolongations of the segments 
@, 4, +--+, a,, and the angles at P;, P,,---, P,,_, right angles. 

> ° ac 4 

Let angle Pil A, = @; then each of the angles P,P,A,, P3P,Ag, - - - 
P,P, An—1 Will also equal 6. 

Moreover let 


PA, = Y= ay, PA, = ", PrA, = eft % P,A,, = Yn: 














1905) PHYSICAL SOLUTIONS OF GENERAL EQUATION OF NT DEGREE 69 


Consider now any one of the similar triangles A,PP,, A\P,P,, - - - 
A,_1P,—1P,, say the triangle A,_,P,_,P, (Figure 5 ). We have 





tang =A _ Pity — Avidyne ~ a 
Py Ayr P,_,A,_; Ye 


or, putting tan 6 = 2, 

Yr = Yr-1 Bt AY. 
Hence, giving to / successively the values 0,1,2, - - -, n, we obtain equations 
[3] of §1, for which we have already shown that 


Yn = Agx™ + ye"? + aye" +--+ Ay ye + Uy. 


In order, therefore, to evaluate this polynomial for a given value of x, say 
x =a, itis only necessary to take 6 = tan —! a, and perform the construction 
indicated. A, will represent the corresponding value of the polynomial. 

If the coefticients a), a), - - -, a, are not all positive, the construction is 
modified in that the angle between consecutive segments a,,, and a,, of the 
line PAyA, -- - A, is taken equal to — 90° whenever a,,, and a, have unlike 
signs, to + 90° when they have like signs. 

The reason for this rule lies in the fact that when we change the sign of 
the middle one of any three segments @,_), “,, @,4,, and the sense of the 
turn which determines its direction, keeping the absolute directions of the 
segments a,_, and a,, , unchanged, the direction of the segment a, with refer- 
ence to each of the segments a,_, and a,,, will have been reversed, so that 
the change of a single sign is accompanied by two changed turns in the line 
PAjA, - ++ Ay. 

Finally the sign of 7’,.A,, = y, is determined by the sign of A,_,A, = a,, 
and y, has the same or opposite sign as a, according as /,,A,, has the same or 
opposite sense as A,,_ ,A,. 

If a is a root of [1], 7A, must equal zero. Now a= tan@. We shall 
therefore know the roots of [1], if we can find the values of @ for which P,, 
coincides with A,,. 

These values of @ can be found in various ways. The simplest way is to 
draw the line PA,.A, --- A} ona piece of tracing paper, prolonging the lines 
AyA;, A;Ag, -- +, A,—1A,, indefinitely in both directions. The tracing is 
then placed over a sheet of millimeter coordinate paper so that P-Ay will coincide 
with one of the ruled lines. This line is taken for the line /P, and the angle 
6 is zero. If now we gradually turn the tracing paper in a positive direction 
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about the point P, increasing the angle @ from zero, it is evident that all posi- 
tions will be reached for which 7, and «1, coincide. After some experience 


these positions can r radily be determined visually. 


3. A dynamic solution. This solution makes use of the principle 
of moments. Suppose we have a set of balances (Figure 6) : 


Syl, Sh, Sy Lhy,*- as Spear, SR, 


with fulerums at the middle points Cy, Ci, C2, ---, Cui» Cn respectively. 
These scale beams are arranged in a vertical plane, successively one above and 


Fia. 6. 





jn 
(la 


RP» 1 














to the left of the other, so that each right hand scale beam, except the lowest, 
may be connected at a distance + 7 from the fulcrum to the right-hand scale 
pan of the beam immediately above it. Let the scale beams be of equal length 
2u. 

Suppose now that weights 


Moy Ay, Ags ° 
are placed in the scale pans at 


Ry, R,, R,, ~ os . = R,, 


respectively, and that 7, denotes the resultant force at any point /?, due to the 
weight a, and the weights a, _,, a,_2, ---, @,, @» in the scale pans of the 
beams below it. 
Then 
Yesr = N+ Opgis [4] 


where .Y is the force resulting at /?,,, from the action of y, at Rt. The dis- 
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tances of #,4,and #, from C;, are v and ~ respectively ; hence by the prin- 
ciple of moments v.Y = y,u, 
or Y= yy F 
a kh . vd 
“ [5] 
where ow. . 
. 
Equation [4] therefore becomes 


Yrer = YE + Opa [6] 


Giving to & successively the values 0, 1, 2, ---,n—2,n— 1, and remember- 
ing that 7%» = a, We obtain once more the equations (3) of §1, from which equa- 
tion [1] follows as before. 

If the weight a, ,, had been placed in the left pan, at S,.,, instead of in 
the right pan at /?,,,, equation [4) would obviously be replaced by 

Yes = X — Oy 41, 

so that for every weight placed in a left pan, the corresponding coefficient in 
equations [3], and consequently also in equation [1], is negative. 

Now since y¥, denotes the resultant force effective at 7?,, the following 
dynamic method of evaluating any polynomial of the form [1], for any real 
value of x, suggests itself. 


We adjust our system of scale beams, x in number, by varying v so that 
u 
v 

+, S, A, weights respectively proportional to the coefficients a,0,, - - -, @, of 
the polynomial, in the right or left pan according as the corresponding 
coefficient is positive or negative. The counterweight required to balance the 
beam S,/?, represents the value of the polynomial for the given value of x. 
y, is positive or negative according as the counterweight must be placed at S,, 
or at 22, in order to balance the system. 

To determine the roots of the equation [2], the scale beams must be 
connected so that the system may be balanced by changing the distance v 
instead of by the use of counterweights. This is accomplished by attaching 
the fulerums C,, C., - - -, C,, to parallel slides, each moved by a separate screw. 
These screws are geared together in such a way that as C, moves through a 
certain distance, C, moves through twice that distance, C,, through x times 
that distance. At the same time each beam slides through the ring by which 
the scale pan at the 7? above it is attached to it. Thus however the distance v 
may be altered, it will have the same value for each beam. 


=a. We then place in one of the pans of each of the beams S,/t,, SR), 
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72 MORITZ . [April 

The pans having been properly weighted, the distance v is gradually 
varied until the system of scales is in equilibrium. For each v that this occurs 
the equation 


furnishes a root of the equation [2]. 

Theoretically v may have any real value negative as well as positive. But 
for practical reasons v is positive and at most equal to unity, hence x is neces- 
sarily positive and not less than unity. But since the roots of 


Agr” a aa-* an*~* as ¢eo & An—1 eo = a, = 0 


are the same as those of [2] with the signs changed, and those of 


— 2? 
2 ar" ° 


a,x" + a,_,2"-'+ 4, fre +art+a=9 [8] 


the reciprocals of the roots of equation [2], the positive roots of [7] with the 
sign changed will give the negative roots of [2] and the reciprocals of the 
roots of [8] whose absolute values are less than unity, will yield those roots of 
[2] whose absolute values are greater than unity. 


4. A solution by hydrodynamics. In the preceding solution the 
balance was brought into equilibrium by changing the lengths of the lever- 
arms at which the weights act. Obviously equilibrium might also be estab- 
lished without changing the distances by a gradual variation of the weights. 

Referring again to the polynomial [1], suppose that at the distances 
My, 4, Any + + +, @,—1 from the fulerum C’ of a scale beam (Figure 7), to the 
left or right of C according as the coefficient is positive or negative, the re- 
spective weights 

Was Weis Waar -- +> W, 


be attached, and at a distance unity from the fulcrum a weight a,, to the left 
or right according as a, is positive or negative. The moment of rotation y,, of 
the system in a positive direction about C’ will then be 


Ay W,, + a Wri + @ Wr—2 Pore oe Gwe W, +4,.- [9] 


Suppose now that each weight W, consists of a vessel containing water 
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and that the shape of the vessel is such that the volume, and consequently the 
weight of its contents, varies as the A*” power of the height to which it is filled. 


FIG. 7. 
Pf ee fj = eee eee ee 
i % ro— Cf — n-l —" 


2 —1— 


6 6 
Ws Wr-2 (res) 


Disregarding the weights of the empty vessels and supposing that they are all 
filled to a common height x we should have 








W,, 


= cx", W,_, = cx"-!,..., Wy = cr, W, = 7, [10] 


where ¢ is an arbitrary factor of proportionality which may be so chosen that 
the height of the water in the vessels is expressed in the same unit as the dis- 
tances from c at which the weights are attached. If the weights are in grams 
and the distances in centimeters, c must be put equal to unity, since the weight 
of a given quantity of water in grams is equal to its volume in cubic centime- 
ters. With ¢ = 1 the values [10], substituted in [9], give us the polyno- 
mial {1}. 

If therefore we could construct vessels having the assumed forms, the 
value of the polynomial [1], for a positive value x of the variable, might be 
found from the counterweight necessary to balance the system after the vessels 
have been filled with water to a common height x. If the units used are grams 
and centimeters, the counterweight at the unit distance will give y, directly, 
and it is positive or ncgative according as the counterweight is required at the 
right or the left of the fulcrum. 

Now it is easy to show that vessels of the assumed form can actually be 
constructed. In fact for every value of / there exists a solid of revolution which 
satisfies the required condition that its volume equals the " power of its 
altitude. 

For supposing that such a solid of revolution exists, take its axis for the 
x-axis, a tangent at its vertex for the y-axis: and let 7 = (x) represent the 
equation of a meridian of the solid. Putting V for the volume cut off from 
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the solid by a plane perpendicular to the x-axis ata distance x, we have on 
the one hand, by the principles of the calculus, 


V= | miypder, 


and on the other hand, by our assumption regarding the form of the vessel, 


V=e'; 


" 
a =| mide: 
J0 


Differentiating this last equation with respect to x we obtain Aa*—! = wy?*, or 


whence 


k—1 


ih 
=W5°*" 1d 
Y \ T U | ) 
Fork =1, 2, 3, 4, 5, we have respectively y = 0.564,(0.798) x4, (0.977) a, 
(1.128) 23, (1.262) 2. 
The meridian sections of the corresponding vessels are shown in Figure 8. 


UL | 


tre re eae PB 
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FIG. &. 











The values of x for which the moment of rotation y¥, is zero, that is, for 
which the system is in equilibrium without the use of counterweights, give 
the roots of the equation [2]. 

These values could be found by gradually increasing the height of the 
water in the vessels, keeping it at the same height in all the vessels, and ob- 
serving the successive heights for which the system balances. The difficulty of 
keeping the height exactly the same in all the vessels is overcome by replacing 
the vessels by heavy solids of revolution of the same shapes and immersing 
these in water to a common depth. By the principle of Archimedes the up- 
ward force of the displaced water is exactly equal to the former downward 
force of the filled vessels. But the direction of the forces having been reversed, 
the weight a, must now be attached at the opposite side of the fulcrum, that 
is at a unit distance to the right of ( if the coefficient a,, is positive, to the left 
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if negative. The gradual immersion of all the bodies to an equal depth is easily 
effected by placing the apparatus in a vessel into which water may be slowly 
admitted. The bodies are adjusted so that their lowest points lie in the same 
horizontal plane. The effect of the weights of the bodies is eliminated by 
carefully balancing the system before attaching the weight a,. Then the weight 
a,, is attached and water admitted to the vessel. The successive depths of im- 
mersion at which equilibrium takes place give the positive roots of equation [2]. 
To tind the negative roots we first transform the equation into another whose 
roots are equal to the roots of the given equation with contrary signs. The 
positive roots of the transformed equation will, of course, give the negative 
roots of the original equation. 


5. An electromagnetic solution, The preceding methods have dealt 
with real roots of equations only ; the present method furnishes all the roots 
of an equation. We shall therefore write the polynomial in the form 


SF (2) = aoz" + ay2z"-' 4 ---a,_,24+ a, [12] 


and seek the values of 
3s + ty 
which will render f(z) = 0. 
Take x + 1 arbitrary constants Aj, Ay, +--+, Ans An41, form the function 
F(z) = (2 — Ay) (z — A) + (2 — An) (Z — Angi), 
and decompose the rational fraction f(z) //'(z) into partial fractions. We get 
. J (2) = M, M,; M,, + M,, +i 


= @& © co o + 


F(z) z—-y | z—-dh —2—ak ZH Aga 





[13] 


where 


My, My, +++) Muy Musi 
are constants given by the relation 


al Fax) 
M, = FOy)? 


F’'(z) denoting as usual the derivative of F(z) with respect to 2. 

We next take some non-conductor of electricity, say a pane of glass or a 
sheet of paper, to represent a plane of complex numbers; on it we select a 
convenient line as an z-axis, and along this line, from a point O as origin, 
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measure the distances OR,, OR, ---, OR,, OR, +, equal respectively to 
Ais Ao + +s Any Angas At the points F,, 2, +--+, Lt, R,, +1 let the horizontal 
plane or pane be traversed at right angles by electric currents whose strengths 
are proportioned to M,, Me, - ++, M,, M4 respectively, with their directions 
downward for positive W's, upward for negative W's. 


The relative current strengths are easily obtained by dividing a single 


circuit into branches by means of wires of uniform thickness whose various 
lengths are inversely proportional to the W's. The strengths of the currents 
carried by these branches vary inversely as the resistances (Ohm’s Law), the 
resistances vary directly as the lengths of the wires, and the lengths of the 
wires vary inversely as the quantities 1; consequently the current-strengths 
vary directly as the quantities .V/. | 

Let 


~ 


z=u+y 
represent any point P in our pane (Figure 9). 


FG. 9. 





By Biot and Savart’s law the force exerted upon a small magnet at P by a 
current of indefinite length and intensity .W, descending perpendicularly at 7?, 
is directly proportional to .W and inversely to p, the distance between 7? and 
P. Its direction is YS, the normal to 27. Hence denoting this foree by F, 
we may write 
_ kM, 
Pp 
where /: is the factor of proportionality. 

Let ¢ denote the angle which the line R/? makes with the z-axis ; then the 
angle which PS makes with /’7, the perpendicular from / to the z-axis, is 
also ¢. The components Y, Y of the form /” in the positive directions of the 
coordinate axes are therefore respectively 


F [14] 


X = F sin ¢, Y = — F cos ¢, 
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whence 
y — hMy y — A M(a — 2) : 
p* ’ , p? ’ [15 ] 


since in the triangle RPT 
: TP y RT «x-—x 
Sin — = * ° ‘Os = = , . 
a at aie es 
Each of the currents at RR, 1, +++, Ry, Ity4, gives rise to x and y com- 
ponents of the form [15]. Hence the resultant force at 7 in the direction of 
the x-axis is given by the sum 
y hMy ks My . [—2 4 May rl May a M, +] ”) 
Pp" Pp Pi P2 Pu Pu +] ’ 


and the resultant in the direction of the y-axis by 


II 


eT AM =) _ yy Me=®) 


3 . —k 
p p 
[Ae —r,) M,(x — Xr) M(t — An) | Mua (@—Ay4)) 
— oo k ; . ee | ; : A, 0 0 0 ke oe —_ + —— ° 
Pi P2 Pr Phy | 
At the neutral points of the electric field we have 
. M _ M(x — 
2 J 0, » 3 cs , at, (16) 
p* p* 


and therefore also 


V(x — } z—r— 
Mie led eS (v A ”) = 0. 
p “ ai 
But 
p* =(xr — r)? + = (x —A+ iy) (x —~A— ty) 


whence 


= (9 Se = *) —> : M ——,. > M = J (2) = () 
p* r—-A+ Ly z—r F(z) 
by [13], and consequently also 
S(2z) = 9. 
It appears, therefore, that the coordinates of every neutral point in our 
plane give a value of z for which f(z) = 0. Moreover for every root of 
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(z) = 0 there exists a neutral point in our field, since every step that led from 
[16] to f(z) = 0 can be reversed. 

The neutral points in question can be determined by various methods. 
The most obvious method is to scatter iron filings over the horizontal pane 
which represents the plane of complex numbers. By slightly tapping the pane 
while the current is flowing, the filings arrange themselves in lines, represent- 
ing the lines of force. The neutral points or points of equilibrium can then be 
easily distinguished. 

If the non-conducting pane is replaced by a metallic plate which receives 
the currents for which .W is positive and discharges them when WV is negative, 


then the neutral points which we have considered are approximately replaced 


(the more accurately the larger the plate) by the nodal points of the system of 
equipotential lines and as such may be located either electrolytically or by 
means of the galvanometer. 


UNIVERSITY OF WASHINGTON, 
SEATTLE, WASH. 














THE BALLISTIC PROBLEM 


By Frank GILMAN 


Tue problem of the motion of a projectile in a medium resisting as any 
power of the velocity was first solved by John Bernouilli in 1718, but he ex- 
pressed the coordinates of any point of the trajectory in a form where a final 
integration is required to determine them numerically ; and, as far as the gen- 
eral problem is concerned, no later mathematician has improved on his results. 

Modern experiments in gunnery, notably those conducted for the British 
government (1865-1880) by General Francis Bashforth, have established the 
fuct that for moderate velocities, ¢. e. from 100 to 900 feet a second, the re- 
sistance of the air to the passage of the projectile is very nearly proportional 
to the square of the velocity : that for greater velocities ¢. ¢. from 900 to 1300 
feet a second it is nearly proportional to the cube of the velocity ; and that for 
high velocities, ¢. ¢. from 1300 to 2780 feet a second it is again very closely 
proportional to the square of the velocity. 

In 1873 General Bashforth published a report of his experiments accom- 
panied with the mathematical theory for the case where the resistance of the 
air is as the cube of the velocity, and an elaborate set of tables based on the 
method of quadratures for the integration of the final equations. 

In 1890 he published a revised report in which the more difficult case 
where the resistance is as the square of the velocity was treated in like manner. 

The same subject has been somewhat differently treated by Niven,* 
Siacci,t and Ingalls. 


In this paper I shall attempt to show that by means of the application of 


the method of Least Squares to the development of functions the differential 
equations for the trajectory when the resistance is as the square of the velocity 


can be integrated without excessive labor, so that they will give results of 


very considerable exactness. 
I shall tirst consider a method ot very great accuracy but requiring much 





* Niven, Proceedings of the Royal Society, London, Vol. 26, 1877. 
t Siacci, Balistique Extérieure; Traduction annotée par P. Laurent; Paris, 1892. 
t Ingalls, Exterior Ballistics, Washington, 1900. 
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labor in its application, and I shall afterwards develop much simpler working 
methods of sufficient accuracy for practical purposes. 
The following* are the well known differential equations of motion for 
the case under consideration : 
dp 


cdx oe emer, (1) 
pyl + p? + log(p+yl+ p*)-C 





: = - (2) 
pyl + p? + log( p + yl + p*) a HE 


edy = sea 


— @p ‘ 
Veq.dt = — —————— : = (3) 


\ C'— pyl + p? — log( p + yl +p?) 


in which 7 denotes the differential of the abscissa and dy the differential of 
dy. , , , 
the ordinate. while p, or a is the tangent of the angle at which the trajec- 
ar 
tory is inclined to the horizontal. 
We have also 
1 2 2 
2ch co? dh, | Poy) + po + log (po + yl + Po) 


' ; q , ; : : : 
in which ¢ = jad being the acceleration of gravity and / the velocity of the 
ve 
projectile when the resistance becomes equal to its weight; 4 = - where vy 
. 29 


is the initial velocity : $o is the initial angle of elevation and py the initial value 
ofp. It is easily seen that C’ is positive and always greater than 


First method. Let us replace the transcendental function pyl + p? 
~ log(p + y1 + p*)which we will represent by /(p) by the rational expression 
a'p + b'p® + ep’ in which the coefficients a’, 6', and c’, are to be calculated by 
the method of Least Squares, the general form of the equation of condition be- 
ing as follows: 


a'p'+ b'p® + ep’ —f(p)=9%. 
Make successively in this equation p = 0.1, p = 0.2, ete. upto p = 1.0, 





* v.'Bartlett’s Analytical Mechanics. 
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thus forming 10 equations of condition. | Forming the normal equations and 
solving them we obtain 


Loga’ =0,.3012124 
Log 4’ = 9.5126277 


Log (—c') = 8.4899010 

A more detailed exposition of this method of developing irrational and 
transcendental expressions is given in an article that was published in the 
Journal of The Franklin Institute for June, 184, entitled "Application of the 
Method of Least Squares to the Development of Functions.” 

The above values of a’, 4', and c’ when substituted in the rational expres- 
sion will make it give values for f(p) whose greatest error will hardly exceed 
one in the fourth place of decimals, while the integrals, cx and cy, are likely 
to be even more exact. 

In formulas developed by the method of Least Squares (especially when the 
number of conditional equations is large) the sum of the positive errors in a 


y 








Go r 
Fis. I 
series of results, will be nearly equal to the sum of the negative errors 
Therefore if the line Oafd represents the course of a given curve with reference 
tothe axes x and y, its location as determined by the method of Least Squares 
is shown by the curve Ocd; and although these curves may deviate sensibly 
from each other at certain points, as at a and c, yet the areas included between 
the lines Oe, ed and each curve will be practically the same. In other words 
the corresponding definite integrals will be nearly equal. 
The differential equations (1) and (2) may now be written as follows: 


dp 


cc'dx _ ' ' vy? (4) 
7 . .. @ ( 
Yo a + 3 3 
a . 5) 
ecl'dy = —— o: ( 
o py + Ps yi e 
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These equations can now be integrated by the method of decomposition 


into partial fractions. 


’ y 


}’ , 
Let / = i m= - , a= sek and observe that / and m are negative and 
c c , 


n positive. 

By Descartes’ Rule of Signs p> + lp} + mp + n = 0 has two positive real 
roots, one negative real root and two imaginary roots; call them a, §, ¥, 
§+ey—1, 6— eV_-l. Then 

p J N Pp 
1 L p+ Q (6) 


p + lp — mp + n p—a of p-B - pP-Y * (p—$)? +e’ 











L p-a 7 1 
i la re ogrs ” "STWTa’ 
Similarly, ; 
a= 5B4+ 312 +m’ 
. 1 
- hy + Bly? + in” 








Clearing ('\) of fractions and writing the coefficients of »‘ and p' equal to zero, 
noting thata —- 8+ y + 26 = 0, we have 


L-M-+-N + P=0,and La + MB + Ny + 278+ YV=0O, 


whence 


P= —(L - M + N), 
Q@=—(1a+ MB + Ny + 2P8%). 
Equations (4) and (5) now integrate into 


, 


c's = Liog(p — a) + Mlog(p — 8) + Nlog(p — 7) +4 4 log[ (p—8)? +] 


+ ie + @ tan—! (* “s °) ‘ (7) 


€ 
= La log(p—a) + MB log(p — B) + Ny log(p — y) 
2P8 ~ ¢ (Ps — Pé _ 
sp ——+ Sbegi(p = 8)? +e] + sah be tan—! (? *) ‘ (8) 
€ 


The constants of integration are determined from the condition that when 
p is equal to its initial value p,, « = 0 and y = 0; but the formulas may be used 
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as they stand without the determination of the constants of integration if pis 
taken between proper limits. 

We will now apply these equations in solving the following numerical ex- 
ample which is given by Professor Bashforth in his treatise on the Chrono- 
graph, page 97. 

Initial angle of elevation ¢, = 30°; initial velocity V = 751 ft. per sec. ; 
diameter of shot = 6.27 in.; weight = 70lbs. Here ¢ = 0.0000339822 
C = 3.455425. _ 

We first find the five roots of the equation 

p+ lp +mp+n=0, 
U 


: : 4 aS a a’ 
in which /=>= — 10.53872354, m =— = —64.7606866, n = — ris +111.84085. 
c c 


The following are the values of the roots : 
a = +1.3779718 
B= +3.6221343 
y = — 4.0293468 


a+B+y 
» 


e = 38 + 1 — y(a + B) — aB = + 5.3255100 


= — 0,4853796 


= — 


Next computing the five constants by the formulas of the preceding page 
we get: 
L = — 0.0093670, M = + 0.0026236, V = + 0.0013514, 
P = + 0.0053920, Q = + 0.01L40848. 
It is evident that the accuracy of the integrals as given by equations 
(7) and (8) will depend on the exactness with which the following equation 


is satistied : 





l 
pyl + Pp + log(p 4 yl+p)y-—C 
1 L M N Pp+Q | 
— a _ ——— + . 9) 
= alse tp at puna * G8 Fe (") 


Substituting in this equation the values of the roots and constants already 
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found, we obtain the following table, which gives the values of the two mem- 
bers of this equation corresponding to ditferent values of p : 


P FIRST MEMBER SECOND MEMBER 


0.05 0.29803 0.29803 
0.10 | 30721 0.30722 
0.15 O1708 0.31704 
0.20 2757 0.32758 
0.25 BBSI5 ; 33896 
BO £85130 35132 
ae 36478 .b6479 
40 | BT95T 37957 
45 39590 B95R9 
£50 | 41404 41403 

| 43436 1.43434 

5727 0.45725 

248334 0.48331 

.51327 0.51325 


4801 0.54801 





In accordance with what has been said in connection with figure 1, it 
follows that the integrals are likely to be even more exact than might be in- 
ferred from the indications of the above table. 

Substituting the values of the roots and constants in equation (7), and 
taking it between the limits p = tan 30° and p = 0, we find x = 6026.8 feet, 
which is the horizontal distance from the origin to where the projectile attains 
its greatest elevation. 

Professor Bashforth finds by quadratures x = 6032 ft., which is 5 feet in 
error. 

Taking equation (%) between the same limits, we obtain y = 1864.4 ft. 
Professor Bashforth tinds y = 1867 ft. In order to tind the range we will as- 
sume with Bashforth that the angle of fall = — (37° 3’), and taking equations 
(7) and (8) between the limits p = tan 30° and p = tan — (37° 3’), we obtain 

v= 11326.1 ft. = range, 

y= — 17.2 ft. 
According to Bashforth, # = 11534 ft., and y = — 6.5 ft., showing an error 
of about 8 feet in the range. 
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The equations (7) and (9) undoubtedly give the coordinates of the tra- 
jectory correct to the nearest foot; but their use involves great labor, and 
they are chiefly valuable as a criterion by which to compare the results of 
simpler methods. 


Second method. If instead of replacing our f(p) = pV 1+ p? 
+ log( p + V1 + p*) by a’p + Ui’ + cp? we assume f(p) = ap + bp®, a and 
» being constants to be determined as before, our method while open to the 
objection that two separate sets of formulas are required one for the ascending 
and the other for the descending portion of the trajectory, is very much less 
laborious in practice than the one we have been considering. 

Since f(— p) = —f(p), if we determine a and + on the assumption that 
p is positive we must use ap — bp? for f( p) when p is negative. 

Integrating the differential equations (1), (2), and (3) on the assump- 
tion that f( p) = ap + bp? we get for the ascending portion of the trajectory : 











1 a+ 2hp — Va? + 460 
ese log + ip 3 = 1 (10) 
eva? + 460 a+ 2hp + va* + 460 
1 . ‘ ax 
y= 3p. 108 (© — ap — bp*) — 97,” (11) 
2h 
t= J sin (sa) ' (12) 
Vbeq Ve + 460 
For the descending portion where f(p) = ap — bp we get: 
r= — 2 — log a - yr — , ifa* > 440, (13) 
cya? — 460 a — 2hp + ya? — 460 
a ( ee =), if a? < 400, (13') 
ey4bC — vio — a 
1 . x ar 
y=- he log(C' — ap + hp”) + ah’ (14) 
op Y b+ 5 ] a 
= “Ee oa 2_ pb + — |. ‘ 


As the constants of integration have not been introduced in the above 
formulas, we must remember in using them to take p between the limits for 
which x, y and ¢ are to be calculated. 
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The following table gives the values of a and 4 corresponding to different 
angles of elevation. 


2 


p t) | log a | log b a a 
0.1 | 5° 43’ | O0.3007161 | 8.6654758 | 1.998555 3.994222 
0.2 | 19 2998279 | &.9559684 1.994472 3.977917 
0.3 42 2983910 | 9.1254675 | 1.987884 3.951682 
0.4 s | 0.2964379 | 9.2440980 1.978964 3.916298 
0.5) 26 34 | 0.2940089 | 9.3342507 | 1.967927 3.872735 
0.6, BO. 58 | 0,2911508 | 9.4060248 | 1.955018 3.822097 
0.7 | 3: — O.2879118 | 9.4649223 | 1.940492 3.765508 
0.8 B& - 2843400 | 9.514827 1.924598 3.704077 
0.9 41 59 | 0.2804821 | 9.5563036 1.907577 3.638850 
1.0 45 © | 0.2763813 | 9.5925171 | 1.889650 3.570776 


! 





In using this table no interpolations will be necessary, as either one of two 
consecutive sets of values of the a and 4 functions can be used without making 
any material difference in the results. The quantities @ and 4 in the table are 
determined in the usual manner by the method of Least Squares, the general 
form of the equation of condition being 

ap + bp? —f(p) =%. 
Make successively in this equation p = 0.01, p = 0.02, ete., up top = 0.10, 
thus forming 10 equations of condition. Solving the resulting normal equations, 
we obtain the log a and log / given in the first horizontal line of the table. 
In the same manner the log a and log 4 found in the second horizontal line of 
the table are formed from 20 equations of condition, and those in the last hor- 
izontal line are formed from 100 equations of condition. 

We will now solve by the above equations and table the numerical exam- 
ple that we have already worked by the first method, and in which C' = 3.455425, 
¢ = 0.0000339822. Since in this example the initial angle of elevation equals 
30°, we take the values of a and 4 from that line in the table corresponding to 
p=.6. Substituting these quantities together with the values of C and c 
in equations (10) (11) and (12) and taking them between the limits corre- 
sponding to p = tan 30° and p = 0, we obtain 

« = 6026.7 ft., 
y = 1864.8 ft., 
t= 10.38 sec. 
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As the angle of fall is 37° 3’, we take for the descending branch the values 
of a and 6 corresponding to p = 0.8. Substituting these in the equations (13’) 
(14) and (15) and taking them between the limits p = tan (37° 3') and p = 0 
we find 


x= 5300.7 ft., 
y = — 1871.1 ft., 
(‘= 11.15 see. 


Range = 6026.7 4+ 5300.7 = 11327.4 ft. 
Time of flight = 10. 38 + 11.13 = 21.51 see. 
If we use the values of @ and 4 corresponding to p = 0.8 for the ascending 


branch of the trajectory, we obtain the following for the x, 7 and ¢ of this 
branch : 


x = 6026. 2 ft., 
Y = 1865. 3 Riss 
C= 10. 38 sec. 


From this it appears that we may generally use the same values of @ and 
} for the ascending as for the descending branch, those quantities being taken 
which correspond to the angle of fall. 

We will illustrate by one more example taken from Capt. Ingalls’ Treatise 
on Exterior Ballistics, page 104. 

A 12-inch service projectile is fired at an angle of departure of 10°, with 
an initial velocity of 1886 ft. per sec.: it is required to find 2), 4, 4, corre- 
sponding to p = 0, and z,, ¥,, 4 for the point where the angle of fall is — 13°, 
both sets of coordinates being referred to the mouth of the gun. 

In this example C = 0.72037, ¢ = 0.0000254985. Take a and 4 from 
that line of the table corresponding p = 0. 2. Substituting in equations (10), 
(11) and (12), and taking p between the limits tan 10° and 0, we find 


a, = 13203.8 ft., 
UF 1294.4 ft., 
t, = 8.46 sec. 


For the descending branch take the @ and 4 corresponding to p = 0.3, 
and substitute in equations (13) (14) and (15), taking p between the limits 
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tan (— 13°) and 0. 
ferred to the vertex of the trajectory : 


Whence 
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We thus obtain the following values of the coordinates re- 


x 


y 
t 


Vy 


Y2 
t. 


Capt. Ingalls working by 


for these quantities : 


v 
" 
X2 
Y2 


= 13203.8 + 9705.3 
1294.4 — 1027.3 
8.46 + 8.32 = 


= 13180.7 ft. 


= 23007.0 ft. 
241.7 ft. 


9705.3 ft., 


— 1027.3 ft., 


8-32 sec. 


1291.8 ft. 


t; = &.44 sec., 
tf, = 16.83 sec. 


9 


, 


, 


— 22909.1 ft., 
= 267.1 ft., 
16.78 see. 


“Niven’s Method” obtains the following values 


If the actual law of the resistance is that of the square of the velocity these re- 
sults of Capt. Ingalls are somewhat erroneous, the x, being 98 feet in error. 


The more exact equations (7) and (8) give the following values for these 


coordinates : 


Third method. 


13205.6 ft., 


1295.5 ft., 


22908.9 ft., 


267.3 ft. 


Equations (1), (2) and (3) integrate at once into 


cL = : log(C — ap), 


re. 
ioe’ 

2 

Vo -t=9 VU — ap, 


C . , 
[p+ | lon (C — ap) ], 


A rougher method which in point of simplicity leaves 
little to be desired, and yet is surprisingly accurate in application is based on 
the assumption that f(p) = PY 1+ p? + log ( p + yl + p*) = ap. 


(16) 


(17) 


(18) 
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and the ascending and descending portions of the trajectory do not need to be 
dealt with separately. 

Ilere again as the constants of integration are not written we must in 
using the equations take the difference between the values corresponding to 
the extreme values of p. 

The value of a@ is to be determined by the method of Least Squares, our 
equation of condition being ap — f(p) = 0; but as we have but one unknown 
f(p) 

p 


The following table gives the values of a for different angles of elevation : 


we have merely to compute the arithmetical mean of the values of 


Pp 1) a p d a 
i 
O05 r 2.00028 O55 28 49 2.03364 
0.10 5 43 2.00124 0.60 30) 58 2.03976 if 
O.15 8 32 2.0272 0.65 33001 2.04635 | ie, 
0.20 1] 19 ? OOLT4 0.70 35 0 2.05339 
0.25 14 2 2.00751 0.75 36 «52 2.06087 | 
0.30 16 42 2.01041 0.80 38 40 2, 06879 
0.35 1 17 2.01408 O85 40 22 2.07711 
0.40 ?1 48 2.OLS1S O90 41 59 2.08585 
O.45 24 14 2.0224 O.95 43 32 2.09497 
O50 26 «34 2.10447 


2.02799 1.00 45 0 

: = t 
Applying this third method to the examples already worked by the first ye 
method and the second method, we get for the range and time in the first ex- 


ample 
x = 11335 ft., 


t= 21.51 sec.: 


and for the coordinates and the time in the second example 


Xe = 22906 * 
266 ft., 


16.78 sec. 
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A THEOREM CONCERNING UNIFORM CONVERGENCE 
By Geo. D. Birk Horr 


Ix this note I aim to prove the converse of the theorem: If the series 
(1) lay (ax)! + | ug(x)] 4 
converges uniformly in the interval ¢ S 2 5 4, then the series 
(2) wy(u) r My(0) cn 
converges uniformly in the same interval, however its terms may be re- 
arranged. * 


The converse theorem is the following : 


Theorem: Jf a series (2) ts uniformly convergent in the interval 


\|A 


a xr <h howerer tts terins may he re-arranyed, then the absolute value series 
(1) is uniformly convergent in the same interval.t 
To show that series (1) is uniformly convergent, is to show that for every 


e there is an .V such that whenever n 2 UV, 
R,(") <€ 


for every point 2 of the interval. 


Here 2, (2) denotes the remainder of the series (1) after the x” term: 
R(x) = es ry Uy 42 (") ee 
and since all the terms of (1) are positive, /?,(2) = 7, (#). 


The proof is indirect, as follows. Suppose that under the given condi- 
tions the series (1) were nof uniformly convergent; then for some fixed e 


* For this theorem, and fora proof that the uniform and absolute convergence of (2) in 
some order does not sutlice to ensure its uniform convergence in every order, see Bécher, 
ANNALS OF MATHEMATICS. ser. 2, vol. 4 (1902-03), pp. 159-160. 

+The question of the validity of this converse theorem was suggested to me by Professor 
E. H. Moore. 

+ This familiar definition of uniform convergence is here repeated at length because it 
is not always given correctly even in the best text-books. Cj., for example, Picard, Traité 
Panalyse, vol. 1, 2d edition, 1901, p. 211; or Goursat, Cours d’analyse, vol. 1, 1902, p. 406. 


(90) 
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UNIFORM CON VERGENCE 


there would be, for every m, some integer n 2 m and some point «,, of the in- 


terval, such that 
R,,(2,,) > €. 


Indeed, since /2,,(x,,)2 ,(“,,), we may say that for every integer m 
there would be some point z,, such that 


RR, (Lm) > €. 





But since (2) is convergent for x = x,,, however its terms may be re-ar- 
ranged, we know that the absolute value series (1) is convergent for xz = z,, ; 
so that we can find an integer r > m such that 2,4 ,(¢m) < Ry(en) — € 
whence /?,,(2,) — Rn+- (tm) > € 


or 
(3) Un +1 (Ly) + Unm+e2 (r,,) reso} Um + r (Ly) > €. 


This inequality may be used as a recurrent formula for the determination 
of an infinite sequence of integers, 


Ny, Wg, My, ++ +, 


as follows: when m= 0, r= m,; then m, +7, = mg: ++ +3 m+ 7) = M44, 
where the r; corresponding to each m; is so determined that (3) holds. 
These integers divide both the given series into sections or blocks, the 7 
section, for example, containing the terms from the (#,; + 1)" to the(a;~,)", 
inclusive. 

Consider a typical section of series (2), namely 
(4) Uma i(lm) + Um ge(Xu) toe e + Buna (2un) 3 
and let S’, be the sum of the absolute values of the terms of one sign in this 
section, and S” the sum of the absolute values of the remaining terms (if 
any). Then from (3), 


mn 
m7 €, 


N ’ 4 S 


m 


so that one or other of the two quantities S}, and S/) must be > «2; we may 


assume that 


(>) Si, > 5) ° 


Now re-arrange the terms of (4) in such a way that all the terms which be- 
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long to S’, come first ; and do the same in each section of the series. The result 
will be a certain re-arrangement of series (2) which may be denoted by 


(6) U(x) + Ui(x) + U3(x) +--+ 


We proceed to show that, contrary to hypothesis, this new series (4) 
cannot be uniformly convergent. 

For, if the first p terms (p S ”) in a typical section of (6) are the terms 
which belong to Sj,, then, by (5), 


= - - € 
On41(2n) * Cn42 (2m) te ep Onsp(Xm) >? 5) ; 


that is, for every .V there are two integers, and m +p(where V<m<m +p), 


and some point z,,, such that 


me 


€ 
» 


- 


, 


, , 
t m(®m) — R m+p (lm ) > 


where /?' denotes the remainder in the series (6); and therefore either 
e , 7 a) >! . 2 ~ 
Ri, (Lm) or else an (Xm, ) Is >e4. 
Hence the supposition with which we started must be false, and the 
theorem is proved. 
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